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A novel attempt has been made to probe the QCD phase boundary by using the experi- 
mental data for transverse momenta of <f> mesons produced in nuclear collisions at AGS, 
SPS and RHIC energies. The data are confronted with simple thermodynamic expecta- 
tions and lattice QCD results. The experimental data indicate a first order phase tran- 
sition, with a mixed phase stretching the energy density between ~ 1 and 3.2 GeV/ fm 3 
corresponding to SPS energies. 
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1. Introduction 



Numerical simulations of Quantum Chromodynamics (QCD) on lattice at finite 
temperature and density predict that at high temperature and/or energy density 
hadrons melt down to a new state of matter, called Quark Gluon Plasma (QGP) 
which might have existe d after a few microseconds of the Big Bang. Lattice QCD 
based calculations ' ' ' ' indicate that the phase transition from a hadronic state 
to a state of deconfined quarks and gluons occurs beyond a critical temperature 
(T c ) and a critical baryonic chemical potential (/i,f ). The order of the QCD phase 
transition and the properties of the QGP state have been a matter of intense debate, 
lately. According to the lattice QCD (1QCD) calculations, the order of the phase 
transition depends on the quark masses and more importantly, on the baryochemical 
potential {hb)- These calculations suggest that the QCD phase diagram has a point 
in temperature-baryonic chemical potential (T — /is) plane where the first order 
phase transition ends. The exact location of this point in terms of T and \ib is not 
known yet, but the 1QCD calculations suggest that it might be within the reach of 
heavy-ion experiments. Since the magnitude of the baryochemical potential at the 
central rapidity region depends on the collision energy (it reduces with the increase 
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in beam energy), by scanning over a broad energy range, it would be possible to 
study the change in the nature of QCD phase transition along the phase boundary 
and probe the critical point. 

A combination of several observables is required in order to study the quark- 
hadron phase boundary and to locate the QCD critical point. In this article, we 
point out that <fi meson may be used as a sensitive probe for studying the phase 
boundary. The <fi meson is the lightest bound state of hidden strangeness (ss), and 
is produced early in the reaction. The <f> meson is least suffered by hadronic re- 
scattering. Relatively long lifetime (~ 45 fm/c) of 4> mesons may ensure that the 
decay products are least affected by re-scattering with other hadrons as the decay 
occurs outside the fireball ™^. On one hand, the theoretical calculations <^" show 
different amount of modification of the <fi meson properties in the medium due to its 
interaction with other mesons in the thermal bath. On the other hand, the results 
extracted from the analysis of the experimental data on the ratio of various hadronic 
species indicate that the (inelastic) interactions of <j> with other hadrons is not 
significant for temperatures below T c . In view of the other systematic effects entering 
to the measurements of the transverse mome ntum distribution of <f>, the change in 
the properties of <fi in the hadronic phase ^-^ may not be significant enough to alter 
the results of the present analysis. The kinematic observables and flow phenomena 
of the <fi mesons are particularly interesting for asserting the equation of state and 
the nature of the phase transition. 



2. Thermal Properties, Bjorken Energy density and 
Deconfinement Transition 

The initial energy density plays a very crucial role in the evolution of the fireball 
which is estimated as follows. The volume of a cylinder of radius R and length dz is 
given by dV = nR dz = ttR 2 t coshy dr\ where t = y/t 2 — Z 2 , is the proper time and 
rj = i rn-rzf, is the space-time rapidity. Identifying r\ with the kinematic rapidity, 
y and considering the volume element of width dy in rapidity space around mid- 
rapidity (y = 0) at the initial time, r one gets the value of the initial volume as 
dV = ttR 2 tq. If the total initial energy within the rapidity window dy is dEx then 
the initial energy density egj, as was first discussed in a seminal paper by Bjorken 
"^ is given by, 

/dE T \ 1 /dN\ 1 

£Bj = \^/^ = W/ <mT> ^' (1) 

where dN/dy is the number of hadrons per unit rapidity and (tot) is the mean 
transverse mass. In the absence of detailed hydro dynamic cal culations, this expres- 
sion has been widely used by various experiments LUiLLLLLalMI to estimate the initial 
energy density. However, the Bjorken energy density given by Eqn. 1 could be much 
smaller than the actual initial energy density produced in nuclear collisions. The es- 
timated energy density in the framework of Bjorken model may be taken as a lower 
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bound I II due to strong reduction, when compared with the initially produced 
energy density from longitudinal hydrodynamic work during the expansion. 

In a first order phase transition scenario, the pressure increases with temperature 
until the transition temperature T c is reached, then it remains constant during the 
mixed phase, and continues to increase after the end of the mixed phase. In a similar 
fashion, the deconfincment transition can be studied by observing the variation 
of average transverse momentum ((px_)) of hadrons as a function of the hadron 
multiplicity at mid-rapidity (dN/dy) ^". The (pr) is expected to reflect the average 
temperature (Tth) an d a flow component which can be related to the initial pressure. 
Similarly the hadronic multiplicity reflects the entropy density of the system. While 
the purely thermal component of (px) cannot be related to the initial temperature 
which remains immeasurable directly above T c , however the flow component in the 
inverse slope can reflect the plateau of the pressure during mixed phase. It has been 
observed that in central heavy ion collisions, the (tot) (= \/"P% + m2 ) of pions, kaons 
and protons as a function of dN^/dy shows a step- like behavior as predicted in^ 
for a wide range of collision energies ' ' \ Hydro dyna mic calculations assuming 
a first order transition could reproduce these data ' ' [ Here, we have made an 
attempt to make the study for 0-mesons which could provide a better understanding 
of the QCD phase boundary, because of the nature of the interaction of </> with 
hadronic matter. 

The tra nsverse mo mentum (pr) spectra of </>-mesons measured for AGS ' ' \ 
SPS 1221231 aiK j RjQc | 24 | 25 | ener gi es have been analyzed in order to obtain the 
inverse slope parameter, T c ff. This effective temperature has contributions from 
both the (random) thermal and the collective motions in the transverse direction. 
It is well known that the inverse slope (T c ff) of the transverse momentum spectra 
of a hadron of mass m can be related to the 'true' freeze-out temperature (Tth) and 
average radial flow velocity (< v T >) at the decoupling surface as: 

T e g = T th + -z m M I ( 2 ) 

The transverse flow extracted from the px spectra of </> meson is particularly useful 
because of its small re-scattering cross-section with hadrons, the </> meson is expected 
to obtain most of its collective flow from the partonic phase which is produced in 
the early stages of heavy-ion collisions. 

Because of long lifetime of 0-meson it can probe the reaction time scales in a 
definite fashion. The px spectra of 0-mesons may reveal quite a few information 
about the evolution of t he sy stem. These can be obtained by making a blast wave 
(BW) model description 1^ of the experimental spectra. For this purpose, we have 
chosen data corresponding to only central collisions of Au-Au or Pb-Pb at different 
energies. Although the BW parametrization of the mx spectra could give a set of 
values for T c ff and < v r >, we have taken the best possible values (with reduced 
X 2 < 1) for our present studies. The results of the BW fit have been shown in Fig.Q] 
and the extracted values of < v r > and T t h are listed in Table I. The table also 
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enlists T c ff and other basic parameters of different colliding systems at different 
collision energies. The similarity of the extracted T t h to the critical temperature 
T c predicted by lattice calculations ^ imply that the <fi meson freezes out near the 
phase boundary and could be used to extract the properties of QCD matter near 
the transition point. Moreover the large values of u r at RHIC energies indicate that 
QGP has undergone substantial radial flow. 

A compilation of measured data for central heavy ion collisions for T c ff, 
for l° w Pt range (0 < Pt < 3.0 GeV/c) and at mid-rapidity has been 
made | 20 | 21 | 22 | 23 | 24 | 25 | 27 | 28 | 2 l as shown in Table L T he resu lts are depicted in 

Fig. [2] as a function of collision energy ^/snn- Both the hadronic and leptonic decay 
channels of (/)-mesons are included in this figure. The T t ^ and v r in Eq. ^ have been 
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Fig. 1. The invariant yield of </> meson as a function of its kinetic energy at mid-rapidity for 
different colliding energies. The top figure shows the results from the STAR collaboration at RHIC. 
The bottom figure shows the results at SPS and AGS energies. The theoretical curves are obtained 
within the framework of blast wave model. 
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Table 1. The table below gives the basic parameters characterizing the colliding systems at different energies, 
viz., centre-of-mass energy, colliding system, centrality, pseudo-rapidity density of charged particles at the central 
rapidity, estimated Bjorken energy density, T e g, radial flow velocity, thermal temperature. The values of v T are 
obtained from blast wave fit. The estimated error on v T is about 10%. 



v /s NN (GeV) 


System 


Centrality 


dN ch 
drj 


e B j(GeV/fm :i ) 


T cB (GeV) 


^r 


T th (GeV) 


Ref. 


4.87 


Au-Au 


0-5% 


120 


1.3 


0.196±0.02±0.02 


0.39 


0.118 


20 21 


6.3 


Pb-Pb 


0-7.2% 


226 


1.0 


0.196±0.02±0.02 


0.43 


0.102 


22 23 


7.6 


Pb-Pb 


0-7.2% 


251 


1.05 


0.237±0.18±0.02 


0.43 


0.143 


22 23 


8.7 


Pb-Pb 


0-7.2% 


272 


1.1 


0.244±0.09±0.06 


0.43 


0.150 


22 23 


12.3 


Pb-Pb 


0-7.2% 


310 


1.25 


0.240±0.08±0.01 


0.43 


0.145 


22 23 


17.3 


Pb-Pb 


0-5% 


426 


2.9 


0.299±0.07±0.01 


0.5 


0.171 


22 23 


62.4 


Au-Au 


0-20% 


472 


3.5 


0.328±0.06±0.02 


0.55 


0.173 


El 


130 


Au-Au 


0-11% 


520 


4.7 


0.379±0.05 ±0.05 


0.63 


0.176 


24 


200 


Au-Au 


0-10% 


670 


1.9 


0.359±0.05±0.02 


0.63 


0.157 


El 


200 


Au-Au 


0-10% 


670 


5.1 


0.376±0.024±0.02 


0.63 


0.174 


[25] 



obtained from the BW description of the transverse momentum spectra of </>. I t is 
necessary to mention here that even at top RHIC energies, it has been observed ^^ 
that a change of impact parameter of the collision doesn't make the required change 
of energy density to observe a change in phase of the produced matter (no step- 
like behavior in < pr > has been seen). So it is imperative to study the variation 
of temperature with energy density or center of mass energy for a fixed collision 
centrality. 

From the experimental measurements of dE^/dy or dN^/dy with (tot), ebj-t 
can be estimated for all centralities, colliding systems and collision energies. In 
the absence of any accurate estimate of the formation time r, we use the quantity 
of eBj.T. Figure [3] shows T% s (top) and T t * (bottom) of (^-mesons as functions of 
ebj-t. In both the figures ([2] and [3]), the inverse slope parameter of the pr-spectra, 
T c g, shows an increase with collision energy/energy density till lower SPS energy is 
reached. A further increase in collision energy leads to an increase of early temper- 
ature and pressure of the system. As a consequence the transverse momenta of the 
produced hadrons and hence the inverse slope parameters increase with collision 
energy. This is followed by the region of constant temperature which corresponds 
to higher SPS energies, where it is expected that the system undergoes a decon- 
finement transition with the creation of a co-existing phase of partons and hadrons, 
which is signaled by a plateau-like structure in the above spectrum. The resulting 
modification of Equation of State (EoS) suppresses the transverse expansion. This 
observation of a plateau-like structure is equivalent to a liquid-gas phase transition 
with involvement of latent heat of the system, characteristic of a first order phase 
transition. Thus the experimental data indicate a first order phase transition, with 
a mixed phase stretching energy density between ~ 1 and 3.2 GeV/fm 3 . At higher 
energies (corresponding to RHIC energies), T^ s further increases with collision en- 
ergy. The EoS at the early stage becomes stiff again leading to increase in early 
pressure and temperature. The thermal component, T t * , also shows a plateau like 
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Fig. 2. Top panel shows the effective temperature of <^-meson {T® s ) and the bottom panel shows 
the 21? as a function of center-of-mass energy (y'SNN) from AGS, SPS to RHIC energies. 



behavior in going toward the RHIC energies. The extraction of T t * involves the flow 
component, v T , which has an ad ded significance. As a function of collision energy, 
v T shows similar behavior as Tjg'^M A similar pattern is suggested for the elliptic 
flow as a function of collision energy I"™. The observation of a mixed phase goes 
in line with the fact that excitation function of various observables show anoma- 
lous behavior or saturation effects starting at lower SPS energies Iaa | j4 | 3o | oo | or | 
This could be related to the onset of deconfinement corresponding to SPS energy 
regime. In addition, the fact that microscopic transport models, based on hadronic 
degrees o f freedom failed to reproduce the observed behavior of kaon inverse slope 
| js | jij | 4in a i g0 indicated the creation of a deconfinement transition corresponding to 
this energy regime. It has also been observed that the inverse slope parameter in 
p + p interactions increases smoothly with collision energy with the absence of any 
plateau- like structure which signify a first order phase transition'^'. An interesting 
review on onset of deconfinement in nucleus-nucleus collisions has recently been 
done by Gazdzicki et ol. ^M 
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Fig. 3. Top panel shows the effective temperature of <^-meson (TzL) and the bottom panel shows 
the 21? as a function of energy density multiplied to formation time (egj-r) from AGS, SPS to 
RHIC energies. 



3. Summary 

In summary, we observe the anomalous energy (density) dependence of transverse 
mass spectra of </>-mesons produced in central heavy ion collisions. The inverse slope 
of the mx spectrum shows an increase in AGS and RHIC energy regime, whereas it 
is almost constant in the intermediate SPS energies. This observation goes in line 
with the previous observations for lighter mesons. This anomaly could be caused by 
the modification of the equation of state in the transition region of the QCD decon- 
finement transition. The present analysis also indicates a first order phase transition 
stretching over an energy density of 1 — 3.2 GeV/fm 3 , which agrees with the 1QCD 
prediction for the critical energy density (e c ~ 0.7 GeV/fm 3 ) for deconfinement. 
This observation is very exciting in view of the observed anomalies in the excitation 
functions of other observables, namely strangeness production, in view of the onset 
of deconfinement in the SPS energy domain. It would have been equally interesting 
to make a similar study on the spectra of Q hyperons, as it is sensitive to the matter 
equation of state at the early stages of collisions. However, at present, enough data 
are not available for the fl hyperons to make a similar study. The ongoing RHIC 
beam energy scan and the same at SPS will certainly help in getting a clear picture 
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to this problem, with the advent of new data sets in the intermediate energy regime. 
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